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І. WSTĘP 


Potencjały sił międzycząsteczkowych, jakie są obecnie brane pod uwa- 
gę, mają wszystkie wspólną cechę: są potencjałami o symetrii kulistej. 
Założenie symetrii kulistej jest przybliżeniem i to dla cząsteczek, które 
takiej symetrii wyraźnie nie wykazują — dość grubym, spowodowanym 
brakiem opracowań zjawisk przepływu dla potencjałów niesymetrycz- 
nych, wobec znacznych trudności zarówno zasadniczych, jak i rachun- 
kowych. | 

Aczkolwiek już przy bardzo uproszczonym potraktowaniu cząsteczek 
jako sztywnych, sprężystych kulek otrzymuje się pewne teoretyczne 
rezultaty, to jednak dopiero potencjały międzycząsteczkowe typu 
Lennarda-Jonesa, czy Buckinghama zastosowane w teorii 
kinetycznej gazów w sposób dużo pełniejszy ujmują dane doświadczalne. 

Przeprowadzone w tej pracy rachunki opierają się na dwóch typach 
potencjalu Lennarda-Jonesa i na zmodyfikowanym potencjale 


Buckinghama: 
8 4 
о = ae (9) 


көк ep : 
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gdzie e = minimum energii potencjalnej, г = odległość dwóch cząste- 
czek, го = odległość cząsteczek, dla której (ro) = 0, rm = odległość czą- 
steczek, dla której Ф(г„,) = £, a” = dodatkowy parametr określający stro- 
mość krzywej potencjału sił odpychających *). 

Zmodyfikowany potencjał Buckinghama zawiera pewną niedogodną 
osobliwość. Oto dla r > 0, (r) > —o0, a zatem przy odpowiednio małym 
Tmax, P(r) osiąga maksimum, a następnie maleje. Analiza jednak wyka- 
zuje, Ze Tmax jest tak małe, że nie wchodzi w rachubę przy zwykłych zde- 
rzeniach cząsteczek. Możemy więc przyjąć słuszność wzoru 3 dla г 2 Tax, 
zaś ф(г) = со dla г < Tmax- 

Stałe €, ro, rm, « są określane osobno dla każdego gazu. Natomiast 
w przypadku zderzeń cząsteczek dwóch składników gazowych I i II, 
przyjmuje się taki sam kształt krzywej potencjału z odpowiednimi sta- 
łymi £12, (To)12, (гъ), %12. 

Stałe го і Tq są zależne. Z przyrównania pochodnej potencjału do zera 
można prosto wyliczyć stosunek tych stałych (tab. 1). 

Tab. 1. Stosunek odległości ro, dla której ф(го) = о (tzn. odległości przy zderzeniu 


o bardzo małej prędkości początkowej) do odległości rm dla której ф(гһ) osiąga 
minimum 


Charakter potencjału 


Lennard Jones (8,4) 0,8409 
sg sił (GR 0,8909 
Buckingham @ = 12 0,8761 
Бе 1 0,8832 
w = 14 0,8891 

т—15 | 0,8942 


II. STAŁE POTENCJAŁU WODORU I DWUTLENKU WĘGLA 


W ostatecznych wyrażeniach na współczynniki gazowe (przewodnictwa 
cieplnego, lepkości, dyfuzji i termodyfuzji), otrzymanych z teorii przy 
założeniu odpowiedniego oddziaływania międzycząsteczkowego wystę- 
рија stałe potencjału €, то, ги, ©’, o których była wyżej mowa. W związku 
z tym zarysowują się раз Кеп postępowania w zależności od wyt- 
kniętego celu. Po pierwsze, mając z doświadczenia określone wartości 

*) Stała termodyfuzji а jest zupełnie inną wielkością, niż parametr а’ wystę- 
pujący w potencjale Buckinghama. 
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danego współczynnika gazowego w danym przedziale temperatur można 
przez porównanie z teorią, tak wyznaczyć stałe potencjału, by zgodność 
teorii z doświadczeniem była jak najlepsza. Po drugie, mając już określo- 
ne stałe potencjału można je wykorzystać przy teoretycznym wyliczaniu 
mierzonych w doświadczeniu współczynników gazowych. 

Celem tej pracy jest określenie stałych €, го, a, przy czym stałe te 
dla wodoru i dwutlenku węgla (oddzielnie) określane będą z zależności 
lepkości tych gazów od temperatury, zaś dla mieszaniny Н, — CO; z za- 
leżności stałej termodyfuzji a od temperatury. Przy wyznaczaniu tych 
stałych w przypadku wodoru i dwutlenku węgla zastosowałem oryginal- 
ną metodę liczenia, natomiast w przypadku mieszaniny tych gazów opar- 
łem się na metodzie В. М. Sriwastawa 1 K. P. Sriwastawa (1). 

Wzór na lepkość jednorodnego gazu według Chapmana (2) jest 
następujący: 

азе 7. 

Lala a 9 (2) 4 
gdzie К = stała Boltzmana, Т = temperatura bezwzględa, 92) = tzw. 
całka zderzeń Chapmana (2). 

Przy opracowaniu poszczególnych przypadków potencjału różni auto- 
rzy wprowadzali swoiste oznaczenia, różniące się między sobą i między 
oznaczeniami wprowadzonymi pierwotnie przez Chapmana. 

Rozpatrzmy kolejno poszczególne przypadki. 


1. Potencjał Lennarda — Jonesa (8,4) 


Dokładna jego forma została wprowadzona przez piszącego w analo- 
gii do tej formy, jaką użyli Hirschfelder, Bird i Spotz (3): 


ABIT Pa 
= BI Е ү. 
о) -() | 
В. С. Jones (4) wprowadza za Enskogiem całkę zderzeń 227”, 
wiążącą się z całką Chapmana О®(2) następującą zależnością: 


1 
2[4kTe] + 
Q2(2) = 3 ro | | Q 2.» 5 


p? 


m: ma m ) 
m, tm Z 
Zatem wzór 4 na lepkość dla jednorodnego gazu po uwzględnieniu 5, 


przyjmie postać 
1 1 
ЕТ\ч /mkT|2 
15 Б T 
Р 


64 ro” Q 7 
Wzoru powyższego nigdzie w literaturze nie spotkałem. 


gdzie u jest masą zredukowaną 


- (dla gazu jednorodnego |» = 
2 
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2. Potencjał Lennarda — Jonesa (12, 6) 


Hirschfelder, Bird i Spotz (3) przy opracowywaniu tego 
przypadku wprowadzili jeszcze inną całkę zderzeń У (п), która wiąże się 
z całką zderzeń Chapmana zależnością i 

a) 2 Ak. 40) 
9 (n) = ro [ү | W (пй) 


p 
Po uwzględnieniu, że u = =. otrzymamy 
0909) = 2r | اد‎ ۳ ИО) 
т 


(SEE 

EE 6 SE AU р 
16 ro? W®(2) 

Wzór ten podają Hirschfelder, Bird i Spotz w swojej pracy. 


а więc 


3. Zmodyfikowany potencjał Buckinghama 


E. Mason (5), który opracował teorię tego przypadku, wprowadził 
całkę zderzeń, którą oznaczył 2*0%(Т*). Jest to po prostu całka zderzeń 
Chapmana podzielona przez wartość tej całki przy założeniu, że cząsteczki 
zachowują się jak sprężyste, sztywne kulki. 

Przy użyciu tej całki zderzeń wzór na lepkość jednorodnego gazu, jak po- 
daje Mason, wyraża się następująco: 

_5 [МВТ js, 1 

ОЖ | ) Мога? 0*22(Т*) 

M = masa cząsteczkowa gazu, К = stała gazowa, № = liczba Avogadro, 


m 
4 


>ы 


Е 


HR 


Wprowadzając stałą Boltzmana, otrzymuje się 


= Н 

SLP. W. ПЕ >. 

10° ra O ONE 

Tabele wartości całki zderzeń podane są u Masona jednak nie dla 
Q*2,20(T*), leczdla 2022 — [т *(1— r) | z Q* 22 (T*) 


a 


1 


Zatem s 
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W tej pracy rozważane będą dwa krańcowe przypadki związane z po- 
tencjałem Buckinghama dla « = 12idlaa = 15 (Mason podaje 
wartości całek zderzeń 2/22 dla a = 12, 13, 14 i 15). 


mk + 
5 т (acz) 


Wówczas ыык ی کے‎ к 9 
3 mk |! 
gi" 

гыл STRZ тт 10 


We wszystkich omawianych wyżej przypadkach potencjału wartości 
całek zderzeń 07), We), Z2 zostały podane nie w zależności od T, lecz 


T 
od Т* = — (R. С. Jones podaje w zależności od !/2(T*) +). Z drugiej 
ЕЗ 
strony we wszystkich wzorach występuje średnica cząsteczek г, (lub г»). 
= 


Powstaje więc problem, jak dobrać najlepiej tę parę parametrów: —i То» 


by zgodność wyliczonych wartości współczynnika lepkości 1 w porówna- 
niu z doświadczenie była najlepsza. 

Jeżeli wykreślić krzywe Тао. = N(T) oraz Nteor. = AN(T*), to okazuje 
się, że są to krzywe monotoniczne rosnące, tak że trudno znaleźć przy- 
najmniej dwa punkty odpowiadające sobie na obu krzywych. Natomiast 


dośw. . 
Te oraz fy(T*) = А-ну" 


(A — pewna stała), obie krzywe wykazują charakterystyczne maksima. 
Jeżeli te maksima zachodzą dla Т, oraz То*, to z ich odpowiedniości na 
obu krzywych otrzymamy 


w przypadku wykreślenia krzywych f(T) = 


к = т, | к: 
Oczywiście, możnaby podać szereg innych działań tego typu, prowa- 
dzących do wyszukiwania punktów odpowiadających sobie na obu krzy- 
wych. Dlatego też otrzymana w ten sposób wartość € była uważana za 
orientacyjną (chociaż w przypadku potencjału (12,6) była zupełnie do- 
kładna). 
Zanim przejdę do omówienia dalszego postępowania, wypiszę dokład- 
ne postaci funkcji f(T*) dla poszczególnych przypadków potencjału. 


*) _ 64го? yjteor. 1 
Przypadek (8,4) pe” ) = ги ГА ЫК ЫЙЫ TRN" 
ИЕ Px 
Przypadek (12,6) { (T*) — „фри а сот 1 
2 


— |2 
— 
09 


[| (жЕ т* wog 
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Przypadek (exp. 6) 


a =12: {т пае: Lesz”; ao эе... 
2 : є\1 /шКк\1 qT* 7, 22 
(5) [> |" «'=12 14 
A = 15: үт» >. 16rm? EE M 1 | 
Те. 
k 2т @’ =15 


Stałe przy funkcjach f(T*) zostały dobrane tak, by te funkcje można 
było określić dla temperatur, w zakresie których podawane są odpowied- 
nie całki zderzeń. Z przytoczonych wzorów 12—15 widać, że kształt tych 
funkcji nie zależy od rodzaju gazu. Odpowiednie wykresy podane są na 


dośw. 
rys. 1—4. Natomiast wykresy funkcji f;(T) = oddzielnie dla Н; i dla 


105}: (T) 


1 {еог. KT 
Rys. 1. Zależność funkcji f SA. EF ا‎ = m dla przypadku potencjału 
2 
Lennarda — Jonesa (12,6) 


CO; zostały umieszczone na rys. rys. 5 i 6. Do tych ostatnich wykresów 
użyto średnich wartości współczynników lepkości z pomiarów dokonanych 
przez różnych autorów (6). 

Na podstawie zależności 11 otrzymano orientacyjne dane stosunków 
dla Н, 1 СО.. 

Wyniki umieszczono w tabeli 2. 


$ 


k 


Na podstawie danych dla danego gazu przy danym potencjale sił 

Е У ym 
międzycząsteczkowych można ze wzorów 6—8 obliczyć го lub rm, jeżeli 
za ņ weźmie się odpowiednie wartości doświadczalne. Wyliczenia takie wy- 
konałem dla różnych temperatur w zakresie od 100°К do 1200°К (śred- 


Wyznaczanie stałych potencjału Н, СО, i mieszaniny Н, — СО, 7 


Rys. 2. Zależność funkcji | (Т9 =А - a= -od T=“ dla przypadku potencjału 
2 
Lennarda — Jonesa (8,4) 


105% (Т") 
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x 
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Rys. 3. Zależność funkcji f (IEA ваз "Tala przypadku potencjału 


Buckinghama przy parametrze а’ = 12 


nio biorąc). Otrzymałem na ogół różniące się w pewnym stopniu warto- 
ści па г, lub ry w zależności tylko od temperatury. Oczywiście, zarówno 
E - =: > P 
Го jak iar mogą się zmieniać wraz z temperaturą. Ale we wszystkich 
rozważaniach prowadzących do określenia całek zderzeń, ze względu na 
i tak skomplikowane rachunki przyjmowano te wielkości za stałe. Dla- 
tego też różnice w otrzymywanych wartościach ro (czy Ty) w zależności 
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1) teor. 
т" 


Buckinghama przy parametrze а’ = 15 


Rys. 4. Zależność funkcji f (TASA od T="T dla przypadku potencjału 
2 


В 
Tab. 2. Wartości parametru ПЕ otrzymane z równości 11 na podstawie wykresów 
1—6 dla wodoru i dwutlenku węgla w przypadku różnych potencjałów 


od temperatury świadczą raczej o mniejszej lub większej przydatności 
danego kształtu potencjału przy zderzeniach cząsteczek danego gazu, niż 
Е 
k 
Powyższe stwierdzenie wziąłem za punkt wyjścia przy wyznaczaniu 


o teoretycznej zależności го 1 od temperatury. 


dokładnych wartości т, i +" Mianowicie takie wartości dla danego gazu 


i przy danym potencjale uważałem za najwłaściwsze, przy których po- 
mocy wyliczone wartości г, jak najmniej zmieniają się wraz z tempe- 
raturą. Średnie wartości г, dla takiego przypadku dawały jednocześnie 
średnicę cząsteczki Н, lub CO). 

= 
k 
potencjału (8,4), oraz dla CO, w przypadku wszystkich rozważanych po- 


W ten sposób zostały określone stałe oraz го dla Н» w przypadku 


tencjałów. Natomiast stałe T dla H} w przypadku potencjału (12,6) 
i (exp-6) określono bezpośrednio z wykresów 2—5, chociaż powody tego 
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Tab. 3. Tabela przedstawia wartości średnicy cząsteczki wodoru wyliczone w róż- 
nych temperaturach w przypadku potencjału Lennarda — Jonesa (126). 
Przy wyliczaniu korzystano ze wzoru 7 oraz wyznaczonego uprzednio parametru 


= ° 
= 37,25 


r (А) dla H, 


nie były jednakowe. Dla potencjału 
(12,6) wartość określona z wykresów 


2—5 dawała nadzwyczaj dobrą nieza- 
leżność го od temperatury (tab. 3), dla- 
tego też spełnione zostały nałożone 


z góry warunki na. W przypadku 
potencjału (exp-6) dla różnych wartości 
Е 
к D . 
temperatury, przy czym im mniejsze 
osiągane różnice A(ro”), tym większy 
był rozrzut procentowy wartości го’. 
Dlatego w tym przypadku poprzestano 
na wartościach otrzymywanych z wy- 
kresów 3—5 i stąd odpowiednio wyz- 

naczano średnią Го. 

Ostateczne wzory otrzymane z za- 
leżności 6—10, na podstawie których wyliczałem re, po podstawieniu 
odpowiednich stałych gazowych, umieszczone są w tab. 4. W przypadku 
potencjału (12,6) brane były dodatkowo poprawki II rzędu (3). 

Na wykresach 7—9 i 11—14 podane są zależności r, od temperatury 
dla wodoru i dwutlenku węgla oddzielnie dla różnych potencjałów. Grub- 
szą linią zaznaczona jest krzywa, na podstawie której wyliczano średnie 
wartości го. Natomiast na wykresach 10 i 15 zebrano uzyskane wyniki. 
Dzięki tej samej skali wykresy te pokazują, w jakich przypadkach poten- 
cjałów wykreślone krzywe spełniają najlepiej nałożony warunek nieza- 
leżności od temperatury. 


otrzymywano różne zależności ro od 


= 


W tab. 5 umieszczono ostateczne wartości xi ry. 


Warto zaznaczyć, że na wykresach zostały przedstawione zależności го, 
a nie го, od temperatury dlatego, że do wzoru na współczynnik lepkości 1 
wchodzi właśnie kwadrat średnicy cząsteczki. Zatem odchylenia krzy- 
wych przedstawionych na wykresach od prostej poziomej wyrażają jed- 
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10x, (Т) 


4 100 200 300 400 500 К 


Rys. 6. Zależność funkcji {= 19 od temperatury dla dwutlenku węgla 
1 T 


nocześnie odchylenia wyliczonych na drodze teoretycznej wartości współ- 

Е 
№ 
wartości doświadczalnych. Dlatego też zgodność z doświadczeniem przy 
obliczaniu współczynnika lepkości jest najlepsza dla takiego przypadku 
Е 
Е 
kres zależności го’ od temperatury jak najmniej odchyla się od linii 
poziomej. 


czynnika lepkości (dla odpowiednio ustalonych i ro) od odpowiednich 


potencjału i przy takich wartościach i ro, dla którego odpowiedni wy- 
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Potencjał 
Lennarda — Jonesa 
(8,4) 


Potencjał 
Lennarda — Jonesa 
(12,6) 


Potencjał 
Buckinghama 
(exp,6) « = 12 
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Tab. 4. 


3 28437" | „at, 1328T* o 
ачар a 
(z) 1 x 10790) 8 1 х 107877 
sl i a A „ 
хе 7<2) (a) аа 7<2) (А) 
1 Х 107%у?!(2) nX 10W (2) 

205,7T mal RA 961,0T Р 

Го =, 7 | r= и мы ПА 


Potencjał 
Buckinghama 
(exp,6) = 15 


8, 225,2 о 1053Т о 
PAY A А 


Rys. 7. Zależność kwadratu średnicy cząsteczki 


+. ең. E 
różnych wartości Е 


1100 °K 


900 


100 


wodoru od temperatury dla 


w przypadku potencjału Lennarda — Jonesa (8,4). 


Zależność wyróżniona na wykresie była podstawą przy wyznaczaniu 


średniej wartości średnicy 
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ro (A) 
= M6 S C= | ЖЕШ А 
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Rys. 8. Zależność kwadratu średnicy cząsteczki wodoru od temperatury dla 
х ą 
różnych wartości x w przypadku potencjału Buckinghama a = 12. Zależ- 


ność wyróżniona na wykresie była podstawą przy wyznaczaniu średniej wartości 
średnicy 


100 300 500 700 300 1100 °K 


Rys. 9. Zależność kwadratu średnicy cząsteczki wodoru od temperatury dla róż- 


= 
nych wartości Е w przypadku potencjału Buckinghama a’ = 15. Zależność 


wyróżniona na wykresie była podstawą przy wyznaczaniu średniej wartości 
średnicy 
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Rys. 10. Zależność ro? od temperatury, z których wyznaczano średnie wartości го 
wodoru w przypadku różnych potencjałów. Użycie tej samej skali daje wyraźną 
miarę porównawczą tych zależności 


Rys. 11. Zależność kwadratu średnicy cząsteczki СО» od temperatury dla różnych 


s 
wartości = w przypadku potencjału Lennarda — Jonesa (8,4). Zależność 


wyróżniona na wykresie była podstawą przy wyznaczaniu średniej wartości 
średnicy 
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Rys. 12. Zależność kwadratu średnicy cząsteczki СО» od temperatury dla różnych 


wartości Ś w przypadku potencjału Lennarda — Jonesa (12,6). Zależność 
k 


wyróżniona na wykresie była podstawą przy wyznaczaniu średniej wartości 
średnicy 


19 
5 1000 1200 °K 


Rys. 13. Zależność kwadratu średnicy cząsteczki СО. od temperatury dla różnych 
= 
wartości Е w przypadku potencjału Buckinghama a = 12. Zależność wy- 


różniona na wykresie była podstawą przy wyznaczaniu średniej wartości średnicy 
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19,5 


15 


Rys. 14. Zależność kwadratu średnicy cząsteczki СО» od temperatury dla różnych 


e 
wartości к w przypadku potencjału Buckinghama a = 15. Zależność wy- 


różniona na wykresie była podstawą przy wyznaczaniu średniej wartości średnicy 


12 200 400 650 в 100 1200 К 


Вуз. 15. Zależności ro* od temperatury, z których wyznaczano średnie wartości го 
dla CO; w przypadku różnych potencjałów. Użycie tej samej skali daje wyraźną 
miarę porównawczą tych zależności 
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e 
Tab. 5. Ostateczne wartości parametrów ro» (w A)i Ж (w stopniach) przy zało- 


żeniu różnych typów oddziaływań międzycząsteczkowych 


(exp,6) 
@ = 16 


37,3 18,9 


220 


Tab. 6. Porównanie danych doświadczalnych z wyliczonymi współczynnika lep- 
kości wodoru w przypadku potencjału Lennarda — Jonesa (12,6). W kolum- 
nie III podane są wyliczenia Hirschfelder, Bird i Spotza przy parame- 


Е 
trach го = 2,968 A, R 33,3, zaś w kolumnie IV wyliczenia dokonane w niniej- 


Е 
szej pracy na podstawie stałych го = 2,926А, E =37,3° 


т°к | qdośw. qteor. HBS 1 teor. (4n)5HBS (41)% 
100 421 416 415 1,2 8 1,5 
200 681 683 688 у 0,3 W 1,0 
300 896 4 896 906 4 0,0 1.1 
400 у 1084 1073 $ 1095 0,3 1,0 
500 С 1255 Г 1237 1266 0,4 0,9 
600 А ` 2 TA 1407 1425 5i 1,0 0,2 

jj 700 1576 1552 1575 1,6 0,0 

۰800 1 720 1 689 ue 1717 Е 1,5 0,2 
900 1 863 1 829 1852 1,8 0,6 

1 000 M 2 002 1958 1982 2,2 1,0 

1100 2137 ý 2 085 2108 2,5 1,3 

1200 К 2205 2205 2 220 0,0 0,7 
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Rozważania powyższe i wyniki otrzymane wskazują, że zarówno 
w przypadku wodoru, jak i dwutlenku węgla najlepszą zgodność z doś- 
wiadczeniem daje potencjał Lennarda — Jonesa (12,6). Niemniej 
na podstawie tego samego potencjału można dobrać w sposób bardziej 
lub mniej właściwy stałe = i ro. Dla porównania przytoczyłem w tabe- 
lach 6 i 7 wartości współczynników lepkości 1 dla Н, i CO, wyliczone 
Е 
® 
potencjału (12,6), oraz dane obliczone przez Hirschfeldera, Bir- 


przeze mnie na podstawie stałych i г, określonych w tej pracy dla 
da i Spotza (3), którzy przyjmowali dla wodoru г, = 2,968 .A, < = 
= 190°. Te teoretycznie obliczone wartości współczynnika 1 porównane 
są z danymi doświadczalnymi. W kolumnie piątej tabel 6 i 7 podałem 
odchylenia procentowe od wartości doświadczalnych, danych, wyliczo- 
Tab. 7. Porównanie danych doświadczalnych z wyliczonymi współczyniikami lep- 


kości dwutlenku węgla w przypadku potencjału Lennarda — Jonesa (12,6). 
W kolumnie III podane są wyliczenia Hirschfelder, Bird i Spotz'a przy 


s 
parametrach ro = 3,996A, ara 190”, zaś w kolumnie IV wyliczenia dokonane w tej 


pracy na podstawie stałych ro = 3,880A, га = 220°. 


TK | qdośw. nteor HSB о 
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nych przez wymienionych wyżej autorów, zaś w kolumnie szóstej tych 
tabel takie same odchylenia procentowe współczynnika lepkości 1 wyli- 
czonego przeze mnie. 

Jak widać z przytoczonych tabel otrzymałem dużo lepszą zgodność 
z danymi doświadczalnymi, szczególnie w przypadku dwutlenku węgla. 


III. STAŁE POTENCJAŁU MIESZANINY H: — СО. 


Parametr — dla mieszaniny Н; — СО; wyliezałem z zależności stałej 
dyfuzji termicznej @ od temperatury. Stała dyfuzji termicznej (która 
bynajmnie stałą nie jest), związana jest ze współczynnikiem dyfuzji ter- 
micznej От i współczynnikiem dyfuzji zwykłej Dı» zależnością 


gdzie x; i x» oznaczają koncentracje dwóch dyfundujących składników 


gazowych. 
Pierwsze przybliżenie stałej a w oznaczeniach Chapmana jest 
a =5(C— 1)g 16 
gdzie 
C= 29 ө (2) 
FAL 17 
591 (1) 


х151 — Xa8> 


E> 2 
х1 Ө, + x2 9, FX1X>Q;2 19 


Sı = МЕ, — 4AM;M; — 3MM: — Mı) 
Qı == E,[6M;* е (5 == 4B)M;?* Ы> ЗМ: МА] 
Qı = 3(М, — М») (5 — 4B) + 4М,М»„А(11 — 4B) + 2Е,Е; 


1) (1) 1) | 
Әт (2). p — 5% (2) — @2(3) 
1 Су 
5912(1) 591» (1) | 19 
2) (2) 
naa (03 p; = № | 
1 маса 1 
5М„® (1) 5М,9 (1) 
= И نک‎ үч 
РКИ M; mı + ma 
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5, oraz Q, powstają z 5, i Qı przez zamianę wskaźników 1 na 2 
i odwrotnie. 

Wskaźniki 1, 2 przy całkach zderzeń © oznaczają, że dana całka obli- 
czana jest dla przypadku zderzeń dwóch różnych cząsteczek mieszaniny 
gazowej, bądź też jednego tylko rodzaju. 

Czynnik g we wzorze 18 prawie nie zależy od temperatury. Przykła- 
dowo, w tab. 8 podałem wyliczone wartości tego czynnika w zakresie 


Tab. 8. Wartości współczynnika g (wzór 18) obliczone przykładowo na podstawie 
potencjału Lennarda — Jonesa (8,4). Tabela pokazuje. że praktycznie 


є 
wielkość g nie zależy od temperatury, co zostało wyzyskane przy obliczeniu —_® 
k 


dla mieszaniny Н, — СО, 


0,531 | 0,527 | 0,525 | 0,527 


temperatur 200°К — 600°К dla potencjału Lennarda— Jonesa 
(8,4). 


Po zlogarytmowaniu wzoru 16 otrzymamy 
loga = log[5(C — 1)] + log(g) 


natomiast logarytmowanie zależności T = T* Ж daje 
log T = logT* + log a 


Jeżeli teraz wykreślimy dwie zależności 


1. fı(logT*) = log[5(C — 1)] 20 
na podstawie danych tablicowych dla danego potencjału, 
oraz 2. f;(logT) = loga 21- 


na podstawie danych doświadczalnych dla Н; — СО», to te dwie krzywe 


powstają jedna z drugiej przez równoległe przesunięcie o log zł wzdłuż 


osi odciętych i o log (g) wzdłuż osi rzędnych. Chodzi teraz o takie dobra- 
nie punktów odpowiedniości, by krzywe najlepiej się nakładały. W tym 
celu wykreślono styczne do obu krzywych pod trzema różnymi kątami 
(rys. 16—20). Punkty styczności na obu krzywych odpowiadające tym 
samym kątom nachylenia stycznych są właśnie szukanymi punktami od- 
powiedniości. Odczytane na osi odciętych i rzędnych przesunięcia tych 


puntków wyznazają log а i log(g). Ponieważ styczne były poprowadzo- 
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ne pod trzema różnymi kątami, przeto ostateczne wartości gt g są 
otrzymane jako średnie z trzech wyliczeń. 

Jeśli chodzi o wykreślenie krzywych teoretycznych (wzór 20) na pod- 
stawie różnych przypadków potencjału, to zadanie rachunkowe w sto- 
sunku do wyliczania współczynnika lepkości było nieco ułatwione. Do 
wielkości C, którą należało wyliczać, wchodzi stosunek dwóch całek (wzór 
17). I chociaż te całki zderzeń mają różne postaci w zależności od zało- 
żonego potencjału, to jednak przy dzieleniu zostaje zredukowany stały 
czynnik, o który się różnią w stosunku do całek zderzeń w oznaczeniu 
Chapmana. Z tego względu określenie wielkości С podane we wzo- 


logá 


Rys. 16. Logarytmiczny wykres zależności 
doświadczalnych wartości stałej dyfuzji 
termicznej a od temperatury (wzór 21). 
Na wykresie wyrysowano trzy styczne 
nachylone pod kątami. których tangensy 
kolejno wynoszą 1/,, 1, 2. Punkty stycz- 
ności były porównywane z odpowiedni- -0,6 
mi punktami krzywych teoretycznych 


Log [5(C-1)] 


0,695 | 0,665 | 0,695 


64,9 


65,3 | 57,4 


Rys. 17. Wykres funkcji (wzór 20) obli 
czanej przy założeniu potencjału Len- 
narda — Jonesa (8,4). Punkty 
styczności w porównaniu z odpowiedni- 
mi punktami krzywej doświadczlej (rys. 


$ . 
16) pozwalają na wyliczenie ы (wyniki 


zamieszczono w tabelce przy wykresie) 
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tog[5(C-1] 


średniaj 


0,662 "m 0,524 | 0,582 


110,1 | 91,4 14,0 91,8 


Rys. 18. Wykres funkcji (wzór 20) obli- 
czanej przy założeniu potencjału Len- 
narda — Jonesa (8,4). Punkty 
styczności w porównaniu z odpowiedni- 
mi punktami krzywej doświadczalnej 


5 
(rys. 16) pozwalają na wyliczenie —? 


(wyniki zamieszczono w tabelce przy 


/kresie) 
Log" wy 
Log [5(С-1]] 
tg a 1/2 | 1 | 2 ARIE 
g 0,740 0,662 0,635 0,679 03 
e 
E 104,5 Ма 76,7 БЕ 


Rys. 19. Wykres furikcji (wzór 20) obli- 
czanej przy założeniu potencjału Buc- -0,4 
kinghama a = 12. Punkty styczności 
w porównaniu z odpowiednimi punkta- 
mi krzywej doświadczalnej (rys. 16) poz- 


e 
walają ną wyliczenie —? (wyniki za- 
k 


mieszczono w tabelce przy wykresie) -05 


06 "0 loq” 
rze 24 nie zmienia się bez względu na to, jakie z omawianych typów 


całek zderzeń doń wstawić. W cytowanej już pracy Hirschfelder, 
Bird i Spotz'a podane są wyliczone wartości wielkości С w zależ- 
ności od temperatury, a Mason podaje tabelę wartości C* = 5с, 

= 
К 
dzonych na podstawie wykresów 16—19 przy różnych potencjałach mię- 
dzycząsteczkowych. 


W tab. 9 podane są wartości g i otrzymane z wyliczej przeprowa- 
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ш średnia! 


g 0,552| 0,582 0,641| 0,592 


| 
e ka 
z [1169 | 948 в (тв 97,1 


Вуз. 20. Wykres funkcji (wzór 20) obli- 
czanej przy założeniu potencjału Buc- 
kinghama a = 15. Punkty styczności 
w porównaniu z odpowiednimi punktami 
krzywej doświadczalnej (rys. 16) pozwala- 


o + + Е + + . 
ją na wyliczenie —? (wyniki zamieszczo- 
k 


no w tabelce przy wykresie) 


= 
Tab. 9. Wartości g i T wyliczone na podstawie wykresów 16—20 dla mieszaniny 
H: — CO; przy różnych potencjałach. Dla porównania w kolumnie czwartej i piątej 


€ 
przytoczono stałe T wodoru i dwutlenku węgla wyznaczone ze współczynnika lepkości 
Potencjał 


(8,4) 


(12,6) 


(ехр,6) œ = 12 


(exp,6) « = 15 


Ostatnia wielkość, jaka pozostała do określenia, to stała гз miesza- 
niny Н;— СО, (w pierwszym przybliżeniu jest to średnia arytmetyczna 
średnic obu cząsteczek). Skorzystamy znowu z doświadczalnie wyznaczo- 
nego współczynnika lepkości (7), tylko tym razem dla mieszaniny Н, — 
СО». Ze względu na znaczną komplikację rachunków ograniczę się do 
wyznaczenia stałej rı w przypadku potencjału Lennarda — Jo- 
nesa (12,6), który, jak wykazały wyżej przeprowadzone rozważania 
i obliczenia, jest najbardziej odpowiedni w rozważanych przypadkach. 
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Dla mieszaniny 50%0—50%/ (stosunki cząsteczkowe) dwóch gazów 
współczynnik lepkości w oznaczenia Chapmana wyraża się nastę- 
pująco: 


z+ = + = А+ +. —2А 
= 3 Mo z 
"miesz. x 10' = 2 2 22 
ŚĆ LR" 
gm stm во бабы Fife 
1x 107 12 x 10% 21). Тә x10° ЗЕ m;max107 
gdzie miesz. — współczynnik lepkości mieszaniny gazowej, 1, — współ- 


czynnik lepkości jednego składnika gazowego (СО?) — współczynnik 
lepkości drugiego składnika gazowego (H>), 


Е: а 
® 
8M;M» 912 (1) 

Porównania z doświadczeniem dokonałem dla temperatury Т = 291°K. 
WÓWCZAS miesz. * 107 = 1453, т • 10° = 1448, п»: 107 = 871 (cgs). Ponie- 
waż do oznaczeń A i E wchodzą całki zderzeń © zależne z kolei od para- 
metru rı», przeto po podstawieniu danych liczbowych do prawej części 
wzoru 22 otrzymuje się wyrażenie zależne od г. Po wykonaniu odpo- 
wiednich rachunków otrzymałem 


ڪڪ 


4, 188 гу, + 31,49 
0,05731 ти, + 2,947 r + 13,58 


W celu ustrzeżenia się od pomyłki liczbowej w rachunkowych obli- 
czeniach ten sam współczynnik lepkości dla mieszaniny dwóch gazów 
wyliczyłem ze wzoru podanego w oznaczeniach Hirschfelder'a, 
Birda i Curtissa (8): 


. 103 (гу А) 23 


= 
Tmiesz. X 10 


7 HEE 
Timiesz. X = X, +Ү, 24 
gdzie 
R, = 2хух+ 2 | 
i 7 UE 1a 
и ° 2 2 
у= m 2x,X, (m, my)” Ч | Хә Ya 
AGGIE 4MM NMa Nə mı 
3 + т 25 
= Ар ja s e) 
Zn 512 ыт ma F ZXyXg od ly $ т 4 x а m: 


Ши у) 
гр» (2) 


м 
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Aż, i О*, są odpowiednimi wyrażeniami Aj i Q; podzielonymi przez 
te wyrażenia dla przypadku cząsteczek uważanych za sprężyste, sztyw- 
ne kulki. 

W oznaczeniach Үү, 2; można dla prostoty rachunków brać zamiast 
mas cząsteczek wyrażonych w gramach odpowiednie masy cząsteczkowe. 
Natomiast w wyrażeniu na т; tak został dobrany czynnik liczbowy, by 
też za m; i mą można było podstawić masy cząsteczkowe dwutlenku wę- 

` а i wodoru. 

Po wstawieniu odpowiednich wartości liczbowych do wzorów 24—25 
otrzymałem 


A 4,268 rj, + 31,51 
0,05727 r, + 2,952 r, + 13,59 


fimiesz. Х 107 10° (гә w А) 26 
Јак widać z porównania ze wzorem 23 zgodność jest najzupełniej za- 
dowalająca, co wyklucza możliwość błędu rachunkowego. 
Do wyliczeń parametru г!» wziąłem zależność na miesz. ze średnimi 
współczynnikami przy rı2: 


4,268 rı + 31,50 
0,05729 гї? + 2,950 ги + 13,58 


Zależność 27 po uporządkowaniu względem potęg г; przedstawia równa- 
nie dwukwadratowe, które ma tylko jeden pierwiastek dodatni. Pier- 
wiastek ten гі = 3,422 A różni się nieco od średniej arytmetycz- 
nej średnic cząsteczek dwutlenku węgla i wodoru wyznaczonych w tej 
pracy. Średnia ta bowiem wynosi 


7) miesz. Х 10* = 


27 


сакто _ 3,880 + 2,926 $ 3.404 А 
гені: 2 


Nie przemawia to oczywiście na niekorzyść przeprowadzonych rachun- 
ków, gdyż taka średnia może być uważana jedynie jako pierwsze przy- 
bliżenie parametru rı2. 


€12 


Parametry ry» i k dla mieszaniny dwóch gazów mogą być jeszcze 


wyznaczone z tak zwanych zasad rekombinacyjnych. Zasady te, przyto- 
czone niżej, mają wszystkie większe lub mniejsze uzasadnienie teore- 
tyczne. Jeżeli formę potencjału sił międzycząsteczkowych przedstawimy 
w ogólnej postaci 


g(r) = Ae т — Cr”, 
to zasady rekombinacyjne mają postać 


sp=/ 4 ° 2 28 
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m= (r +m) 29 
н-на) u 
Ар =УА: · Аз 31 
Св =V С.С? 32 


W szczególności, zasady 30, 31 i 32 stosował Mason, zaś zasady 28, 30 
і 32 В. М. Sriwastawa і К. P. Sriwastawa (1) dla potencjału 
(exp-6). W przypadku potencjału Lennarda — Jonesa (12,6) pier- 
wsze dwie zasady pozostają bez zmian, zaś pozostałe przyjmą postać 


ACZ й 
ано Zły 

1 
Сиз ry = РГА 34 
© |2 гу = гут) ИЕ, €3 ` 35 


Stosując do tych zasad dane, otrzymane oddzielnie dla wodoru i od- 
dzielnie dla dwutlenku węgla, można nietrudno otrzymać parametry 212 
i гү; mieszaniny Н, — СО». 

W tab. 10 umieściłem w ten sposób otrzymane wartości parametrów 
8, ; s Өү. F А 3 
кїт? przy zastosowaniu różnych wzorów rekombinacyjnych. 

Z tabeli wynika sprawa trochę nieoczekiwana, że najlepszą zgodność 
z doświadczalnie wyznaczonymi parametrami 22 i г dają najprostsze 


zasady rekombinacyjne 28 i 29. 


e 
Tab. 10. Wartości parametrów =. i го mieszaniny Н: — CO; otrzymane z różnych 


wzorów rekombinacyjnych przy potencjale Lennarda — Jonesa 42,6) 


Zasady rekombinacyjne Dane 


z doświad- 


(28),(29) | (29),(35) (33),(35) (34),(35) chała 
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ZAKOŃCZENIE 


Wyniki przeprowadzonych obliczej wskazują, że do interpretacji zależ- 
ności od temperatury współczynników lepkości wodoru i dwutlenku wę- 
gla najlepiej nadaje się potencjał Lennarda — Jonesa (12,6). 
Różnice jakie otrzymuje się między wyliczonymi na podstawie tego po- 
tencjału współczynnikami (przy odpowiednio dobranych parametrach r» 
i fa ), a danymi doświadczalnymi dla Н, i СО, wyrażają się zaledwie dzie- 


siątymi częściami procentu i to w dużym zakresie temperatur. 


Natomiast jest rzeczą charakterystyczną, że mający najlepsze uza- 
sadnienie teoretyczne potencjał Buckinghama nie nadaje się do 
interpretacji danych doświadczalnych współczynników lepkości H і СО». 


Stałe r, i E jakie się otrzymuje na podstawie tego potencjału znacznie 


odbiegają od tych, jakie się zwykle przypisuje cząsteczkom Н, і СО». 

W przypadku mieszaniny H, — СО» przy interpretacji danych ter- 
modyfuzji na podstawie rozważanych przypadków potencjału, otrzyma- 
e 
k 


= 
nimi różnica (szczególnie między wartością + wyliczoną na podstawie 


ne wartości są wszystkie do przyjęcia. Aczkolwiek występuje między 


potencjału Lennarda — Jonesa (8,4), a pozostałymi wartościa- 
mi), to jednak jest ona taka, że da się uzasadnić odmiennym charakte- 
rem rozważanych potencjałów. 

Wreszcie należy odnotować, że przy wyliczaniu stałych potencjału 


rı i "айа mieszaniny Н, — СО, ze stałych określonych osobno dla Н, 
i СО», najlepsze wyniki dają najprostsze zasady rekombinacyjne: 


паа rı Fr: 2 с сч 
B= TZW OMB 812 = | Е, 2 
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PE3IOME 


В настоящей работе подвергнуто анализу применение интермоле- 
кулярных потенциалов типа Леннарда — Джонсона (12,6) 
и (8,4), а также Бекингема для а = 12 иа = 15 к интерпретации 
зависимости коэффициентов липкости водорода, двуокиси углерода 
и смеси Н — СО», равно как постоянной термической диффузии этой 
смеси от температуры. Применение интегралов столкновений в случае 
потенциала Леннарда — Джонса (8,4) к исчислению коэффи- 
циента липкости дало автору возможность найти неизвестную доселе 
формулу относительно этого коэффициента в рассматриваемом случае 
потенциала. Констатировано, что наибольшее согласование с опытом 
в этих пределах дает потенциал Леннарда — Джонса (12,6). 


€ 
Для этого потенциала определены постоянные J» A также го в случае 


Н,, СО», равно как смеси Н; — СО.. Эти постоянные отличаются He- 
сколько от постоянных, полученных Гиршфельдером, Бер- 
дом и Споцом. 

Установлено очень хорошее согласование теоретических данных 
липкости полученных на основании потенциала Леннарда — 
Джонса относительно постоянных определенных в этой работе, 
с экспериментальными данными. 


SUMMARY 


The applicability of the intermolecular potentials of Lennard — 
Jones (8,4) and (12,6) and the modified Buckingham types for 
о/ = 12 and a =15, has been investigated to explain the dependence on 
temperature of the viscosity and thermal diffusion constant of the H2 — 
CO; mixture and of the viscosity of pure Н» and СО). A new formula for 
the coefficient of viscosity has been derived with the help of the collision 
integrals of the Lennard-Jones (8,4) potential. 

The results of the experiment corresponded most closely to the 
Lennard-Jones (12,6) potential. The values of its potential para- 


meters > and r, were found to be somewhat different from those 
calculated by Hirschfelder, Bird and Spotz. 
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